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the  mapping  A{z).  Both  functions  are,  in  general,  nonlipschitzian.  Thus  nonstandard 
techniques  are  required  in  the  analysis  of  their  variational  properties.  In  order  to  give 
some  insight  into  the  problem,  let  7;  (f  k-»  (C"  be  analytic  and  consider  the  eigenvalues 
of  A(Q  =  ^(7(C))  near  zero.  It  is  well  known  that  these  eigenvalues  are  given  by  so 
called  Puiseux  series,  that  is,  series  in  fractional  powers  of  (^  with  the  smallest  such 
power  being  greater  than  or  equal  to  -.  Our  results  are  a  consequence  of  certain 
properties  of  these  series.  The  derivation  of  these  properties  depends  on  a  classical 
computational  scheme  due  to  Newton  [11]  known  as  the  Puiseux-Newton  diagram. 

Our  method  of  analysis  requires  the  introduction  of  a  notion  of  directional  differ- 
entiability that  depends  on  analyticity  in  the  following  way:  For  w:  C"  •— ►  R,  define 
w''(z;):C''  >-^  RU{±oo}  by 

w  {z:d)=      mf     lim  mf ,  (6) 

76r(z,<i)     do  c 

where 

r(z,d)=  {r- €'-*€' \y  is  analytic  ,  7(0)  =  2,  and  7'(0)  =  d}.  (4) 

The  superscript  h  in  (3)  is  used  to  emphasize  that  w'^(z;d)  depends  only  on  the  holo- 
morphic  curves  in  C  that  pass  through  z.  This  notion  of  directional  differentiability 
shares  many  of  the  properties  of  other  such  notions  that  have  recently  been  developed 
in  the  literature  on  nonsmooth  analysis  [5,9,10,14].  Understanding  the  relationship 
between  (3)  and  these  other  notions  is  important  for  the  development  of  a  calculus. 
However,  we  defer  the  discussion  of  these  issues  to  a  future  work.  At  present  we  con- 
centrate on  the  evaluation  of  (3)  for  the  functions  a  and  p. 

Let  'Wci'^]  be  the  set  of  monic  polynomials  in  A  whose  coefficients  are  analytic 
mappings  from  G  C  €"  to  C  If  G  =  C",  we  simply  write  'W[A].  Observe  that  E(i) 
is  by  definition  the  set  of  roots  of  a  polynomial  in  'W[A].  For  this  reason  we  begin  in 
Section  2  by  studying  the  properties  of  a  and  p  when  it  is  assumed  that 

E(z)  =  {A6C:P(2,A)  =  0},  (5) 

where  P  £  'W[A].  This  approach  has  its  limitations  since  in  general  we  only  obtain  a 
lower  bound  on  the  value  of  the  directional  derivative  (3).  Nonetheless,  this  lower  bound 
can  be  shown  to  be  sharp  in  certain  nondegenerate  situations.  Let  Ao(z)  and  TZo{z) 
denote  the  elements  of  E(z)  that  achieve  the  maximum  value  in  (1)  and  (2),  respectively. 
Then,  in  the  context  of  matrices,  the  polynomial  results  are  most  meaningful  when  all  of 
the  eigenvalues  in  either  ^^0(2)  or  TZoiz)  are  nonderogatory,  that  is,  their  multiplicity  in 
both  the  characteristic  and  minimal  polynomial  for  A{z)  coincide.  For  a  nonderogatory 
eigenvalue,  the  characteristic  polynomial  contains  all  of  the  essential  information  about 
the  eigenvalue.  This  is  the  case  in  which  our  lower  bound  is  most  likely  to  be  sharp.  The 
opposite  extreme  is  when  all  of  the  eigenvalues  in  either  Ao{z)  or  7vo(2)  are  semisimple 
(nondefective),  that  is,  when  the  multiplicity  of  these  eigenvalues  in  the  characteristic 
polynomial  is  greater  than  or  equal  to  one  and  yet  their  multiplicity  in  the  minimal 
polynomial  is  precisely  one.    In  this  case,  the  characteristic  polynomial  contains  the 


least  amount  of  information  about  the  eigenvalues.  It  is  in  this  latter  case  that  the 
lower  bounds  that  we  obtain  for  the  directional  derivatives  are  most  likely  not  sharp. 
On  the  other  hand,  one  can  precisely  evaluate  the  ordinary  directional  derivative  in 
the  semisimple  case.  This  is  done  in  [12].  Thus  the  cases  for  which  precise  results  are 
unknown  are  those  where  at  least  one  of  the  eigenvalues  in  either  Ao{z)  or  TZoi^)  is 
both  derogatory  and  defective.  Nonetheless,  the  lower  bounds  that  we  establish  do 
provide  a  great  deal  of  information  about  when  these  directional  derivatives  attain  the 
value  +00. 

In  the  discussion  that  follows  certain  statements  are  sensitive  to  the  domain  of  the 
variable  being  discussed.  Thus,  in  order  to  avoid  confusion,  we  obey  the  following 
convention  concerning  the  labeling  of  the  variables  (^,  e,  and  2:  C  will  always  represent 
a  complex  scalar,  e  a  real  scalar,  and  z  a  vector  in  C". 

2      Roots  of  Polynomials  with  Complex  Coefficients 

Let  P  €  W[A].  We  begin  by  studying  the  differential  properties  of  the  function  g  at 
2°  €  C"  .  From  [2,  pp.  376-381]  there  exists  a  neighborhood  G  of  z°  in  €"  on  which 
P  has  the  unique  representation 

P=     U    ^'>"  (6) 

Ak€S(.!») 

where  /it  G  'Hci'^]  for  each  At  G  ^{z°)-  The  polynomials  /i*  can  be  taken  to  have  the 
form 

/it(A,  z)  =  (A  -  A*)"  +  Ckiiz){X  -  Ai)"-i  +  •  ••  +  Ci<»(2),  (7) 

where 

ct:iz°)  =  0  for  j  =  l,...,h 

and  tk  is  the  multiplicity  of  the  root  At.  We  begin  with  the  case  in  which  P  has  only 
a  single  root  Aq,  with  multiplicity  to,  i.e. 

P(A,  2)  =  (A  -  Ao)'»  +  ci(z)(A  -  Ao)'"-^  +  •  ■  •  +  c,„(2),  (8) 

where 

Cj(2°)  =  0  for  j=l,...,to  . 

We  will  return  to  the  general  case  at  the  end  of  this  section. 

Our  first  objective  is  to  understand  the  behavior  of  g  along  analytic  curves  in  C 
passing  through  2°.  Let  7:  (C  1— »  C"  be  an  analytic  curve  satisfying  7(0)  =  2°.  Compose 
each  Cj  with  7  to  obtain 

^(A,  7(C))  =  (A  -  Ao)'"  +  /?i(C)(A  -  Ao)'''-^  +  ■  ■  •  +  /?,„(C)  =  0,  (9) 

a  polynomial  equation  in  A  with  analytic  coefficients  /?,(C)  =  Cji7(0),  satisfying 

0j{O)  =  O,j=l,...,to. 


Figure  1 


We  may  write 


m)  =  pf\-^p?^('+ 


where,  for  example, 

/?{')  =  c',(z'')7'(0).  (10) 

In  the  discussion  which  follows  we  restrict  (  to  a  nontrivial  real  interval  [0,  to]  and  write 
f  in  place  of  C  to  emphasize  this  restriction. 

As  has  already  been  noted,  it  is  well  known  (e.g.  [2],  [7])  that  the  roots  of  (9)  are 
described  by  series  in  fractional  powers  oft.  These  series  are  commonly  called  Puiseux 
series,  since  it  was  Puiseux  [13]  who  established  their  convergence;  however  they  were 
derived  formally  by  Newton  two  centuries  earlier  (also  see  [15,  Chapter  1,  Section  2] 
and  [11,  page  88]  for  examples  and  applications).  We  obtain  the  results  we  need  by 
making  use  of  a  diagram  devised  by  Newton  for  the  purpose  of  calculating  coefficients 
of  Puiseux  series. 

Let  $j  =  p\  ''  be  the  first  nonzero  value  in  the  sequence  {/JJ  ,/?j  ,...}■  By  defi- 
nition, ij  >  l,j  =  l,...,<o-  If  0}(c)  is  identically  zero,  take  £j  =  oo;  also,  since  the 
coefficient  of  (A  -  Aq)'"  in  P{X,  e)  is  one,  take  ^o  =  0,  $0  =  1-  Now  plot  the  values  ij 
versus  j,  and  consider  the  lower  boundary  of  the  convex  hull  of  the  points  plotted.  Let 
Sj  be  the  slope  of  the  line  on  [j,j  +  1]  forming  part  of  this  boundary,  j  =  0, . . . ,  <o  —  1- 
Clearly  l/<o  ^  sq  <  si  <  ■  ■  ■  <  St„-\.  Figure  1  shows  the  diagram  for  the  following 
example  (taken  from  [2]): 

to  =  3;  Ao  =  0;  /?i(e)  =  c;  /?2(f )  = -(  -  (^\  /?3(0  =  f '  +  Sf^ 


We  have  £0  =  0,  £1  =  1,  £2  =  1,  £3 


2,  and  so  sq  =  si  =  ^,  S2  =  1- 


Now  consider  the  following  "Ansatz"  argument.  Suppose  a  root  of  (9)  is  to  have 
the  form 

A(f)-Ao  =  aeP  +  ---  (11) 

where  a  is  nonzero  and  p  is  the  smallest  power  of  e  in  the  expansion  for  this  root. 
Substituting  (11)  into  (9),  we  need 

(a'of'op  ^  . . .)     ^     (^je'.  4....)(a'o-i^(io-i)P^...)+... 

+     (/Jto-if''"-'  +  •  •  •)(«''  +  •••)  +  (/?<„f''°  +  ••■)  =  0. 

The  terms  involving  the  smallest  powers  of  e  are  among  the  terms 

a'Of'-'P,  /?,a'''-'£''+"°-i)P,...,/?,„_ia£''°-'+P,  /?,„f''o.  (12) 

For  cancellation  to  take  place,  at  least  two  terms  with  the  same  smallest  power  of  e  must 
appear.  Equivalently,  p  must  equal  one  or  more  of  the  slopes  sq,  . . . ,  St^-i  defined  by 
the  Puiseux-Newton  diagram.  The  following  discussion  will  apply  to  a  particular  choice 
of  such  p.  Define  b  and  y  by  p  =  sj  =  . . .  =  Sb+g-i,  so  that  the  line  in  the  diagram 
with  slope  p  passes  from  the  point  (b,£i,)  to  the  point  (6  +  g,£i,^g).  Cancellation  of 
the  coefficients  of  the  terms  with  the  smallest  powers  of  c  in  (12)  requires  a  to  be  the 
root  of  a  polynomial  equation  with  degree  g,  with  leading  term  0i,a^  and  constant  term 
0b+g,  and  with  an  additional  intermediate  nonzero  term  for  each  point  {j,£j)  lying  on 
the  line  in  the  diagram  with  slope  p,  where  b  <  j  <  b  +  g.  Now  let  p  =  q/f,  where 
q,f  are  relatively  prime  integers.  By  definition,  p  is  an  integral  multiple  of  l/g,  so  g 
is  an  integral  multiple  of  /,  say  g  =  mf.  It  is  then  clear  from  the  diagram  that  of 
the  g  —  I  abscissa  values  j  between  6  and  b  +  g,  only  every  /th  value  is  a  candidate 
for  the  intersection  of  the  line  with  a  point  with  integer  coordinates.  Consequently  the 
polynomial  of  degree  y  in  a  reduces  to  a  polynomial  of  degree  m  in  a-' ,  which  we  may 
denote  by  Q{r).  The  conclusion  is  that  the  given  value  of  p  is  associated  with  g  roots 
with  an  expansion  of  the  form  (11),  with  a  taking  the  values 

r;^/V,    h=l,...,m,j^l,...f  (13) 

where  the  r^  are  the  m  roots  of  Q{r)  =  0,  r^'  is  the  principal  /th  root  of  r^  and  w  is 
the  principal  /th  root  of  unity. 

Completing  the  example  given  above,  we  see  that  the  two  values  for  p  are  sq  =  «i  = 
5  and  S2  =  1-  In  the  case  p  =  ^  we  have  fc  =  0,^  =  2,  /  =  2,  m  =  1,  with  Q{r)  =  r—  1, 
so  the  possible  values  for  a  are  ±1,  giving  the  Puiseux  series 

A(f)- Ao  =  ±f^  +  ■•■. 

In  the  case  p  =  1  we  have  b  =  2,g  =  1,/=  l,m=  1,  with  Q{r)  =  r—  1,  so  the  only 
possible  value  for  a  is  1,  giving  the  Puiseux  series 

A(e)- Ao  =£+•••. 


The  subsequent  terms  in  the  series  may  also  be  calculated  by  repeating  the  process. 

The  following  result  is  the  key  to  our  analysis.  It  is  a  generalization  of  [3,  Theorem  1]; 
see  also  [8].  The  lemma  provides  information  about  the  coefficients  of  (9)  when  it  is 
assumed  that  roots  of  (9)  lie  in  the  half  plane  Re  yo(r  — Aq)  <  0  infinitely  often  as  e  |  0. 

Lemma  1  Consider  the  polynomtal  equation  (9),  with  roots  given  by  one  or  more 
Puiseux  series  of  the  form  (11).  Let  yo  £  C  and  suppose  that  there  exists  eo  >  0  such 
that  all  the  roots  X{i)  of  (9)  satisfy 

Re  yo(A(e)  -  Aq)  <  6f  +  o{e)  (14) 

for  some  sequence  {f*}  with  t*^  j  0.   Then 

Reyo/?r'>-<o<5,  (15) 

Rey2/?(i)>o,  Imy2/j(i)  =  o,  (16) 

/?j"  =  0,    i  =  3,...,<o.  (17) 

Here  (16)  is  understood  to  be  vacuous  if  to  =  1. 

Proof  The  coefficient  /3i(f)  is  the  sum  of  the  differences  Aq  —  A(f)  over  the  roots 
A(f)  of  (9);  thus  (15)  follows  from  (14),  letting  e  — ►  0.  The  other  results  follow  from 
the  Puiseux-Newton  diagram  as  follows.  Consider  the  Puiseux  series  corresponding  to 
p  =  soi  the  smallest  possible  value.  In  order  for  (14)  to  hold  either 

(i)p>  1  (e.g.,  if/ =1),  or 

(ii)  p  =  i    /  =  2,  and  Re  yQr^  =  0  for  /i  =  1, ...,  m,  where  the  r^  are  the  m  roots  of 
Q{r),  with  Q{r)  taking  the  form 

Q(r)  =  r'"+4^V'"-'  +  ...  +  4::;).  (18) 

No  other  cases  having  p  <  1  are  possible  due  to  the  splitting  of  the  roots  as  described 
in  (13).  In  both  cases  (i)  and  (ii),  p  >  5,  so  ^j  —  0  for  j  =  3,  ...,<o  from  the  Puiseux- 
Newton  diagram.  In  the  case  p  >  1,  we  also  have  y^j  —  0-  I"  the  case  f  —  2,  observe 
that  the  condition  Re  yQr^  =  0  is  equivalent  to  the  two  conditions  Re  y^ri,  <  0  and 
Im  ylr^  =  0.  Now,  since  — /?2     is  the  sum  of  the  roots  of  Q{r),  (16)  follows.  q 

We  apply  this  result  to  the  evaluation  of  q''(z°;  d)  where  a  is  given  by  (1)  with  E(i) 
and  P  defined  in  (5)  and  (8),  respectively.  First  observe  that  we  can  replace  the  limit 
infimum  in  (3)  by  limit  since  the  perturbed  roots  of  the  polynomial  (9)  are  given  by 
Puiseux  series  of  the  form  (11)  for  some  non-negative  rational  number  p.  Therefore, 
this  limit  always  exists  and  can  only  take  the  value  +00  if  it  is  not  finite.  Consequently, 
o''(z°;  •):  C  I-*  R  U  {+00}  and  is  given  by 

a'^{z;d)=      inf     Hm^Mlllz.^.  (19) 

We  now  state  the  main  result  of  this  section. 


Theorem  2  Let  P  be  given  by  (8)  and  choose  d  £  C".  If  any  one  of  the  conditions 


is  violated,  then 
otherwise 


Re  c'2(z°)(i  >  0,  Im  4(2°)^  =  0, 
cr(z°)d  =  0,    j  =  Z,...,to 


a''(z°;d)  =  +oo; 


Q''{z°;d)>--Rec\iz°)d. 
to 


(20) 
(21) 


(22) 


Moreover,  if  the  rank  of  c'{z°)  is  to,  where  c.C  >-►  (C'"  is  given  by 


c{z)  = 


ci{z) 


then  equality  holds  in  (22)  whenever  (20)  and  (21)  are  satisfied. 
Proof    Suppose  +oo  >  6  >  a''(r°;  d).  Then  there  is  a  7  €  T{z°,d)  such  that 

a(7(c))  -  a(.-°) 


lim  ■ 


<S, 


(23) 


or  equivalently, 


a(7(€))-a(z'')<6f  foree[0,eo], 


q(1) 


for  some  cq  >  0.  Let  /?]''  =  c'j{z^)d  for  j  =  1, 2,  ••  •  ,<o  as  in  (10).  By  invoking  Lemma 
1  with  j/o  =  1.  we  see  that  (20)  and  (21)  must  be  satisfied.  Thus,  if  any  one  of  these 
conditions  is  violated,  we  must  have  a'^{z°;d)  =  +00.  By  letting  6  I  a'^{z°;d),  we  also 
obtain  from  Lemma  1  the  inequality 


O) 


a''iz";d)>--Re0\ 

'0 

Let  us  now  suppose  that  the  rank  of  c'(z'')  is  to-  We  need  to  establish  equality  in 
(22).  Clearly,  we  need  only  consider  the  case  in  which  (20)  and  (21)  hold  for  j/o  =  1- 
In  this  case,  equality  follows  if  we  can  exhibit  a  curve  7  G  r{z°,d)  such  that 


^^^a(7(e))-a(.°)^_]_^      (,) 

no  (  to       '■ 

Consider  the  coefficients  of  the  powers  of  (A  —  Aq)  in  the  polynomial 


(24) 


((A  -  Ao)  +  ^0">-\X  -  Ao  +  i(/?^^'c)'  +  ^C)(A  -  Ao  -  i{/3i'\)i  +  ^Q 
10  to  10 


=     (A  -  Ao)'°  +  0['\(X  -  Ao)'"-'  +  (4"C  +  C»(C'))(A  -  Ao)'°-'  +  . . .  ,    (25) 
=     (A  -  Ao)'°  +  t',(C)(A  -  Xo)'"-'  +  t;2(C)(A  -  Ao)'''-^  +  . . .  , 

where  i  =  \/—\.  Note  that  these  coefficients  satisfy  (16)  and  (17)  with  j/o  =  1-  Also 
note  that  if  7  can  be  chosen  from  (4)  so  that  (9)  has  these  coefficients,  then  (24)  is 
satisfied  for  this  curve  and  the  proof  is  complete.  We  now  show  that  this  can  indeed 
be  done. 

Define  FiC'+'n-f"  by 

F{z,0  =  c{z)-v(0, 

where  v.C  t— >  C'°  is  the  curve  whose  component  functions  are  the  coefficients  of  the 
powers  of  (A  —  Aq)  in  the  polynomial  (25): 

HO  =  [MO,  f2(C), .  •  • ,  MOf  =  <:c'(z'>)y'{0)  +  O(C'),  (26) 

where  the  second  equality  follows  from  the  definition  of  v.  Let  /  C  {l,...,t/}  be 
such  that  the  matrix  c'j(z°)  is  nonsingular  and  set  J  =  {I, . .  .,u}\  I.  By  the  implicit 
function  theorem  [6],  there  is  a  neighborhood  G  C  C"'""*"'  of  ((z'')j,0)  and  an  analytic 
mapping  7:  G  1— ►  C'°  such  that 

F{jizj,0,zj, 0  =  0, 

for  all  {zj,0  €  G  with 

j{{z°)j,0)  =  {z%. 

Furthermore, 

7'((z°).,0)  =  -Mz'r'[c'jiz%-v'm  (27) 

Define  7:  C  >-+  C"  by 

'rjiO  =  (z°)j+Cdj,  (28) 

and 

7/(0  =  7(7^(0,0-  (29) 

Then  for  all  ^  sufficiently  small  c{f{0)  =  '^(0-  Consequently,  with  this  choice  of  7, 
the  polynomial  (9)  is  precisely  the  polynomial  (25).  Moreover,  from  (26)-(29)  we  have 

7^(0)  =  dj 

and 


7K0)  =  7'(7J(0),0) 


dj 
1 


=  di, 


so  that  7'(0)  =  d.  Thus,  7  G  r{z°,d).  This  concludes  the  proof  of  the  theorem.  g 

Let  us  now  return  to  the  general  CEise  and  recall  the  factorization  (6).  Our  results 
in  this  case  follow  directly  from  the  special  case  (8). 


Theorem  3  Let  a  be  given  by  (1)  with  E  and  P  as  defined  in  (5)  and  (6),  respectively. 
Choose  deCr  If  for  some  Xk  €  Ao{z°)  =  {At  G  E{z°)  :  Re  At  =  a{z°)}  any  one  of 
the  conditions 

Re  42(2°)rf  >  0,  Im  c[2{z°)d  =  0,  (30) 

ci^.(z°)<i  =  0,  j  =  3,...,<t,  (31) 


is  violated,  then 
otherwise, 


aHz°;d)  =  +oo;  (32) 


aV;'^)>max{-5^i^£!l^:AtGA(^'')}.  (33) 

'it 

Moreover,  if  the  vectors  {c'i^j(z°)  :  Xk  €  Ai{z°,d),  j  -  \,...,tk}  are  linearly  indepen- 
dent, where 

then  equality  holds  in  (33). 

Proof  The  proof  is  almost  identical  to  that  of  Theorem  2.  The  primary  difference  is 
that  now  all  of  the  roots  in  Ao{z°)  contribute  to  the  value  of  a*'{z°\d).  In  order  to  see 
this  observe  that  the  inequality  (23)  implies  that 

Re  (A(7(e))  -  At)  <  be  for  c  €  [0,  fo]  and  A*  €  Ao{z°), 

for  some  €o  >  0.  Consequently,  (33)  again  follows  from  Theorem  1. 

In  order  to  establish  equality  in  (33),  it  is  again  sufficient  to  exhibit  the  existence 
of  a  curve  -^  ^T(z° ,d)  such  that 

^.^a(,(e))-a(z)^  ^^^_Rec^  _  ^^  ^  ^^^^o^j  (34) 

We  generate  such  a  curve  precisely  as  in  the  proof  of  Theorem  2  except  now  we  must 
choose  the  curve  from  V{z°,d)  so  as  to  match  the  coefficients  in  (25)  for  each  At  € 
A\{z°,d).  Just  as  before,  it  is  the  linear  independence  of  the  gradients  {c'ki{z°)  :  At  £ 
^i(z°,d)}  which  guarantees  that  this  can  be  done  via  the  implicit  function  theorem. 
Moreover,  it  is  clear  that  we  need  only  match  the  coefficients  for  At  G  A\{z°,d)  since 
these  are  the  dominant  first  order  terms.  D 

3      Eigenvalues  of  Complex  Matrices 

Let  W[(C"',  C"^"]  denote  the  set  of  mappings  from  €"  io  €"^"  each  of  whose  compo- 
nents is  an  analytic  map  from  (C"  to  <C.  Let  A  ^HiC ,€"""*].  In  this  section  we  again 


study  the  differential  properties  of  the  function  o  given  in  (1)  with  the  multifunction 
T,(z)  defined  to  be  the  spectrum  of  A(z),  that  is,  T,(z)  is  given  by  (5)  where 


P(X,  z)  =  det[XI  -  A(z)] 


(35) 


is  the  characteristic  polynomial  for  A{z).  One  can  now  apply  the  results  of  the  previous 
section  to  obtain  differential  information  about  a  since  P  G  'W[A].  However,  by  itself, 
this  result  is  not  completely  satisfactory  since  it  does  not  describe  the  relationship 
between  A  and  the  terms  c'i.Az'^)d  appearing  in  (33).  In  this  section  we  describe  this 
relationship,  making  use  of  results  from  [3]  which  in  turn  depend  on  work  of  Arnold 
[1].  Since  our  description  depends  on  the  Jordan  decomposition  of  A^°^  =  A(z°),  we 
need  to  introduce  the  notation  necessary  for  this  discussion. 

Suppose  yl'°'  is  a  matrix  with  eigenvalues  Aj, . . . ,  A^,  having  (algebraic)  multiplici- 
ties <!,...,<,,,  respectively.  Let  the  Jordan  form  of  A^°^  be  given  by 


ylW  =  575-1 


vhere 


J  = 


Jk  = 


Ji 


and  the  Jordan  block 


Jki  = 


Jkx 


A*     1 


J, 


Jknik 


1 

Ai 


has  dimension  n^/.  We  have 

"ii  + l-nimt=<i,      k=\,...,r). 

If  mic  —  1,  At  is  said  to  be  a  nonderogatory  eigenvalue,  while  if  mt  =  <i,  i.e. 
nti  =  •  •  ■  =  rifcmi,  —  1>  ^k  is  said  to  be  semisimple  (nondefective). 

Definition  4  Define  the  jth  generalized  trace  of  a  square  matrix  A,  denoted  by 

tr'j'yl, 

as  the  sum  of  the  elements  on  the  diagonal  of  A  which  ts  j  —I  positions  below  the  main 
diagonal.  Thus  one  obtains  the  ordinary  trace  in  the  case  j  =  I  and  the  bottom  left 
element  of  the  matrix  in  the  case  that  j  is  the  dimension  of  the  matrix.  If  j  exceeds 
the  dimension  of  A,  take  tr'-'M  =  0. 
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Theorem  5  Let  A  e  ■W[(C^  €"'*"]  and  choose  2°  £  C .  Suppose  that  A^"'  =  A{z°) 
has  Jordan  form  as  described  above.  Define 

For  each  q  =  1, . .  .,t/  partition  S~^  A\  S  conformally  with  the  partition  of  J  and  denote 
its  diagonal  block  corresponding  to  Jk  by  Bqk,  k  =  l,...,ri,  with  each  Bqk  having 
diagonal  blocks  Bqki  corresponding  to  Jki,  /  =  1,  ■  -  • ,  "it-   Then,  for  k  =  1, . .  .,t), 


cUz°)  = 


,  forj  =  l,...,tk, 


(36) 


where  the  functions  Ckj  are  as  given  in  (7)  for  the  factorization  (6)  of  P(X,z)  = 
det[XI  -A(z)]  at  z  =  z°. 

Proof  Let  <f  e  €"  and  7  e  r(z°,d)  and  define  A  =  Aoj.  Given  M  e  €"><"  denote 
by  Mki  that  block  of  M  which  conforms  to  the  block  Jki  of  J.  By  [3,  Theorem  4]  we 
have 


c'kji^°)d    =     £tr(-')(5-'i'(0)5)„ 

mj, 

=     ^tr(^>(5-'^'(0)d5)t, 
(=1 

mi,  V 

=     Y.'r'^\S-'{Y,d,A^^^)S)ki 


/=1 


«=1 


q  =  \         (=1 

for  j  =  !,...,<*.  Since  this  holds  for  all  d  eC,  the  result  follows. 


D 


Thus,  given  the  Jordan  form  of  A{z°)  together  with  A'iz")  it  is  possible  to  evaluate 
a  lower  bound  for  a''{z°;d)  by  combining  Theorems  3  and  5.  We  do  this  in  the  following 
corollary. 

Corollary  6  Let  the  assumptions  of  Theorem  5  hold.    If  for  some  Xk  G  Aq{z°)  any 
one  of  the  conditions 

V  mi  "  "In 

Re  ^  d,  ^  tr  (2'fl,i,  >  0,  Im  ^  d,  ^  tr  (^^B,ki  =  0, 

5=1         1=1  «  =  1         '=1 
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5=1 


;=i 


IS  violated,  then 
otherwise, 

Q''{z°;d)>  max 


{ 


a''{z°;d)  = +00; 
ReE:=i'f,E™\trn)B,„ 


A*  G  Ao(z°) 


(37) 


(38) 


Moreover,  if  the  vectors 


for  Xk  e  Ai(z°,d)  and  j  =  1, . . .  ,tk, 


(39) 


are  linearly  independent,  then  equality  holds  in  (38). 


It  should  be  observed  that  if  any  eigenvalue  Ai  G  A\{z°)  is  derogatory,  then  the 
vectors  in  (39)  cannot  be  linearly  independent.  In  order  to  see  this  note  that,  for  at 
least  one  j  between  1  and  <* ,  j  exceeds  the  dimension  of  all  of  the  blocks  J^i  making  up 
Jk,  and  hence  the  corresponding  vector  is  zero.  Thus,  if  ^1(2",  d)  contains  a  derogatory 
eigenvalue,  then  the  sufficiency  condition  of  Theorem  3  is  not  satisfied.  In  this  case, 
one  should  not  expect  to  obtain  equality  in  (33).  Indeed,  in  the  case  where  Aq(z°) 
contains  only  semisimple  eigenvalues,  the  resolvent  theory  for  eigenvalue  perturbations 
yields  the  following  result  [7],  [12,  Lemma  3.5]. 

Theorem  7  If  Ao{z°)  contains  only  semisimple  eigenvalues  of  A{z°),  then  for  every 
d€  C  the  ordinary  directional  derivative  a'{z^\d)  exists  and  satisfies 


a'(z°;d)  =  a''(z°;d)-      max        max  Re  Al,  , 
where  A'j.,,  /  =  1, . .  .,tk  are  the  eigenvalues  ofY^''^  d^Bqk- 


(40) 


Therefore,  if  r;  =  1  and  Ai  is  semisimple,  then  the  right  hand  side  of  (38)  is  the 
average  of  the  real  parts  of  the  eigenvalues  of  Eo=i  '^qBqi,  whereas,  a^{z^;d)  is  the 
maximum  of  these  quantities. 

4      The  Spectral  Radius 

Let  us  now  shift  our  attention  to  the  study  of  the  function  p  defined  in  (2).  We  study 
the  differential  properties  of  p  in  the  same  manner  sis  we  studied  these  properties  for  a. 
That  is,  we  first  consider  the  case  when  the  spectrum  near  z°  G  €"  is  given  as  the  set 
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of  roots  of  a  polynomial  of  the  form  (8)  and  then  extend  this  result  to  the  general  case. 

To  this  end,  consider  the  directional  derivative  p''(z°;d)  defined  by  (3).  As  was  the  case 

for  the  directional  derivative  a''(x;d),  the  limit  infimum  in  the  definition  of  p''{z°;d) 

can  be  replaced  by  limit.  This  is  justified  in  the  same  way  as  it  was  for  a^{z°;d),  that 

is,  by  considering  the  splitting  behavior  of  the  eigenvalues  under  perturbation.  Thus, 

we  may  write 

h,     j\            •   f      I-      P(7(0)  -P(^)  (AU 

p  (z:d)=       mf      lim .  (,41J 

76r(^»,ii)  eiO  £ 

where  r{z°,d)  is  defined  in  (4),  and  p''(;°;  •):  C"  ^^  RU{+oo}.    Continuing  as  in 
Section  2,  we  begin  with  the  following  key  result  for  the  case  in  which  (8)  holds. 

Theorem  8  Lei  P  be  given  by  (8)  near  z°  G  €"  and  choose  deC.    We  will  consider 
two  cases:  p(z°)  =  0  and  p{z°)  ^  0. 

1.  Assume  that  p(z°)  =  0  so  that  Aq  =  0.  In  this  case  we  have 

\   +00  ,  otherwise. 

2.  Assume  thai  p{z°)  ^  0  and  consider  the  conditions 

Re  \oC2{z°)d  >  0,  Im  XoC2(z°)d  =  0,  (43) 

c'j(z°)d  =  0,    j  =  3,...,<o.  (44) 

Then 

o'iz'-  d)>{   ^^  tl^'^^'"^'^'  -^'  ^oc\{z')d\     ,  «/  (43)  and  (U)  hold,      ^^^^ 
P  (     1    ;  _  1^   ^^^  ^  otherwise. 

Here  it  is  understood  that  the  function  C2  is  identically  zero  if  to  =  1. 
Moreover,  if  the  rank  of  c'{z°)  is  <o,  where  €:€"  i—  C*°  is  given  by 

'  c,(z)  " 

c{z)  =      ; 

.  c,(z)  J 

then  equality  holds  in  (42)  or  (45)  depending  on  whether  p(z°)  =  0  or  p{z°)  ^  0  holds, 
respectively. 

Proof    In  this  proof  we  will  continue  to  use  the  notation  of  Section  2.  Let  7  6  r(z  ,  d), 

set  /?^''  =  c^(2°)7'(0)  =  c'j{z°)d  for  j  =  1,2,    -^^o  as  in  (10),  and  let  A(e)  be  one  of 
the  roots  (11)  of  (9).  Since 

|A(e)|2  -  |AoP=  2Re  Ao(A(£)  -  Ao)+  |A(0  -  Ao|^  (46) 
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a  necessary  condition  for 

Sc  +  o{e)>\X(e)\-\Xol  (47) 

or  equivalently, 

26|Ao|e  +  o{€)>|A(Ol'-|AoP,  (48) 

is  that 

6|Ao|€  +  o(€)>ReAo(A(f)-Ao).  (49) 

It  follows  from  Lemma  1  that  if  either  (16)  or  (17)  with  j/o  =  Aq,  or  equivalently,  (43) 
or  (44),  do  not  hold,  then  inequality  (48)  cannot  hold  for  any  (5  £  R.  Since  this  is 
independent  of  7  6  r{z^,d),  p''(z°,d)  =  +00  if  any  one  of  (43)  or  (44)  are  violated 
regardless  of  the  value  of  ^(2"). 

Let  us  now  suppose  that  (43)  and  (44)  hold,  i.e.  (16)  and  (17)  hold  with  j/o  =  Aq. 
Then  for  every  6  >  p''(z°,d)  there  is  a  7  €  r{z°,d)  such  that 

■.n.'"*»-'"-°'<..  (50) 

no  e 

or  equivalently, 

pijie))-p(z°)<6eioT(e[0,co], 

for  some  to  >  0.  Therefore,  inequalities  (48)  and  (49)  are  satisfied. 

As  was  observed  in  Section  2,  the  roots  of  the  equation  (9)  are  necessarily  Puiseux 
series  of  the  form  (11).  Lemma  1  and  inequality  (49)  imply  that  either 

(i)  the  exponent  p  is  greater  than  or  equal  to  1  corresponding  to  series  of  the  form 

A(e)  =  Ao  +  a€  +  o(£), 

with  a  E^  C  possibly  taking  the  value  zero,  or 
(ii)  the  exponent  p  equals  ~  corresponding  to  m  pairs  of  roots  of  the  form 

A(f)  =  Ao  ±  (rhf)^  +  (a±)f  +  o{£), 
where,  as  with  a  in  (i),  a±  takes  different  values  corresponding  to  different  roots. 
By  summing  over  all  the  roots,  we  find  that 

-^l'U^a  +  ^(a+  +  a_).  (51) 

Moreover,  substituting  the  expressions  for  A(e)  given  in  (i)  and  (ii)  into  (48)  yields 

26  |Ao|  t  >  (2Re  Aoa)e  +  o{e)  (52) 

in  Ccise  (i)  and 

26  |Ao|  €  >  (|r,|  +2Re  Aoa±)e  +  0(c)  (53) 
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in  case  (ii).  Now  summing  these  inequalities  over  all  the  roots  gives  the  inequality 

m  

26to  |Ao|  c>[2j2  H\  +2Re  \o(J2  «  +  II(^+  +  ^-))]^  +  "(')'  (^"'^ 


h=l 


where  the  sums  without  explicit  indexing  indicates  summing  over  all  the  roots,  while 
the  factor  2  appearing  in  the  front  of  ^"=i  |rh|  reflects  the  fact  that  there  are  two 
roots  for  each  h=  1, 2,  ■  ■  • ,  m.  Now,  as  explained  in  Lemma  1, 

m  m 

h=l  h=l 

Combining  this  inequality  with  (51)  and  (54)  gives 

6\Xo\>U\Pi'^\-^eXo0['^)+'-^.  (56) 

to  f 

Letting  £  J.  0,  we  obtain  the  inequality 

*|Ao|>^(l/?^'^|-ReAoM'^  (57) 

to 

We  now  consider  the  two  cases  p{z°)  =  0  and  p{z°)  ^  0  separately.  If  p(z°)  =  0, 
then  inequality  (57)  implies  that  /3^''  =  0.  Thus,  p''{z°;d)  =  +oo  unless  ^^^  =  0  for 
j  =  2,  ■  ■  • ,  <o-  Furthermore,  observe  that 

for  all  e  e  [0,eo]  where  the  sum  is  taken  over  all  branches  A(€).  Letting  6  [  p''(z°;d) 
yields  (42).  On  the  other  hand,  if  p{z°)  ^  0,  then  (45)  follows  immediately  from  (57) 
by  letting  (5  ip''(2°;d). 

Next,  suppose  that  the  rank  of  c'{z°)  is  to  and  that  (43)  and  (44)  hold  (otherwise 
equality  in  either  (42)  or  (45)  is  trivially  satisfied).  We  will  only  consider  the  case 
p{z°)  i^  0  since  the  case  p(z°)  =  0  follows  in  a  similar  manner.  Moreover,  in  the  case 
to  =  1,  we  take  /?^^'  =  0  as  usual.  Consider  the  coefficients  of  the  powers  of  (A  -  Aq)  in 
the  polynomial 


((A 


_  Ao)  +  aO'°-'[(A  -  Ao)  +  /^^A+  rf][(A  -  Aq)  -  y^>A+  ^^1 


=  (A  _  Ao)'"  +  /?i^'C(A  -  Ao)'°-^  +  (/?^'^C  +  0(C'))(A  -  Ao)'"-^  +  .  •  •  ,  (58) 

where  a  and  r  are  defined  by  the  expressions 

.=  li[|/?(')|-ReAoM^)]|^ 
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and 


T=\[0['^-ito-2)a]. 


These  coefficients  are  chosen  so  that  not  only  are  (16)  and  (17)  satisfied  with  j/o  =  -^Oi 
but  also  the  expansion  of  each  root  satisfies 

The  proof  now  follows  the  argument  given  in  the  proof  of  Theorem  2,  that  is,  one  uses 
the  rank  condition  and  the  implicit  function  theorem  to  establish  the  existence  of  a 
curve  7  G  r{z°,d)  such  that  (9)  has  the  same  coefficients  as  (58).  With  this  choice  of 
7  equality  holds  in  (45).  g 

The  main  theorem  of  this  section  now  follows  easily  from  Theorem  8.  It  is  derived 
from  Theorem  8  in  the  same  way  that  Theorem  3  was  derived  from  Theorem  2  and  so 
its  proof  is  omitted. 

Theorem  9  Let  P  have  the  representation  (6)  where  each  /ijt  ts  given  by  (T),  choose 
d  6  C",  and  define  no(z°)  =  {Xt  €  T.{z°)  :\Xk\=  p{z°))-  We  consider  the  two  cases 
p(z°)  =  0  and  p{z'^)  ^  0  separately. 

1.  Suppose  p{z'^)  —  Q  so  that  TIq{z°)  =  T,{z°)  =  {Aj}  andi\  =  n.  If  any  one  of  the 
conditions 

c[jiz°)d  =  0,j  =  2,...,n, 

is  violated,  then 

/(z°;d)  =  +oo;  (59) 

otherwise, 

p\z°;d)>^\c\,{z°)d\.  (60) 

2.  Suppose  p{z°)  ^  0.  If  for  some  Aq  6  TZo{z°)  any  one  of  the  conditions 

Re  Xlc't2{z°)d  >  0,  Im  xlc't2i^°)d  =  0,  (61) 

4j(z°)d  =  0,  j  =  3,...,<i,  (62) 

p''(2°;d)  =  +oo;  (63) 

p''{z°;d)  >  max{*(At)  :  X,  eno(z°)},  (64) 


is  violated,  then 

otherwise, 

where 


tkP(2") 


with  cic-2  understood  to  he  the  zero  map  if  it  =  1- 
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Moreover,  if  the  vectors  {c[j(z°)  :  A/t  G  7^l(^°,<f),  ;  =  l,...,<t}  are  linearly  indepen- 
dent, where  7^l(^^<i)  is  the  set  of  Xk  €  7^o(^°)  such  that  ^{Xk)  =  max{*(0  :  ^  G 
TZo{z°)}  when  p{z°)  ^  0  and  TZi{z°,d)  =  Y,{z°)  otherwise,  then  equality  holds  in  (60) 
or  (64)  depending  on  whether  p{z°)  =  0  or  p(z°)  ^  0  holds,  respectively. 

For  the  case  in  which  P  is  given  by  (35)  one  can  apply  Theorem  9  in  conjunction 
with  Theorem  5  to  directly  obtain  a  result  in  terms  of  matrices.  But,  we  again  caution 
that  if  any  of  the  eigenvalues  Xk  G  TZi{z°,d)  is  derogatory,  then  the  linear  independence 
condition  used  to  establish  equality  in  (60)  and  (64)  is  not  satisfied.  Thus,  in  this  case, 
one  should  not  expect  equality  to  hold.  Indeed,  when  Tlo(z°)  contains  only  semisimple 
eigenvalues,  we  obtain  the  following  result  as  an  immediate  consequence  of  [12,  Lemma 
3.4]. 

Theorem  10  Suppose  P  is  given  by  (35).  If'Ro(z°)  contains  only  semisimple  eigen- 
values of  A{z°),  then  for  every  d  E  C  the  ordinary  directional  derivative  p'{zo;d) 
exists  and  satisfies 

{maxtg-RoC^")  maxi<i<„j  |A't,|  «/p(^o)  =  0, 

_ 
^j^maxtgTCoCz")  maxi<,<„^  Re  A^A;,       tf  p{zo)  ^  0. 

(65) 
where  A;,,  /  =  1,. .  .,U  are  the  eigenvalues  o/ 51^=1  d^Bgk  and  the  matrices  B,k  are 
defined  in  Theorem  5. 

Therefore,  if,  for  example,  p(z°)  ^  0,  7^o(^°)  contains  only  the  single  element  Ai 
and  Ai  is  semisimple,  then  the  right  hand  side  of  (64)  is  the  average  of  the  values 
jijRe  AiA'j,,  l=l,...,tk,  whereas  p''{z°;d)  is  the  maximum. 
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